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Abstract 



o. 

2 ' We use a recently found parametrization of the solutions of the inverse Frobenius- 

CJ ■ Perron problem within the class of complete unimodal maps to develop a Monte-Carlo 

approach for the construction of one-dimensional chaotic dynamical laws with given 

^ ' statistical properties, i.e. invariant density and autocorrelation function. A variety of 

j^ ■ different examples are presented to demonstrate the power of our method. 



1 Introduction 

Chaotic dynamical systems show a rich diversity of possible behaviour with respect to their 
statistical properties. In recent years much work has been devoted to the understanding of 
how these statistical properties emerge from the dynamics. The most popular area for such 
investigations are unimodal 1-d maps [|T], 0, |^, Q . They became widespread due to two main 
advantages: their dynamics can be calculated efficiently enough to make extensive numerical 
investigations and the variety of statistical properties within this class of dynamical laws is 
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very large. Certain features of special systems can even be investigated analytically. The 
interplay between the dynamical and statistical behaviour becomes more transparent if we 
consider the inverse problem, i.e. given the statistical behaviour of a system how could 
we extract relevant information on the possible dynamics. This subject has recently been 
addressed and discussed in the literature P, H, ^ |], |, |^ and from a communications point of 
view in refs. |Tl|, . The purpose is then to construct an one-dimensional fully chaotic (and 
ergodic) dynamical law for a given invariant density or/and time autocorrelation function. 
Given the latter we can on the one hand calculate all expectation values, i.e. statistical 
averages, of measurable observables depending exclusively on the single dynamical variable. 
On the other hand the given autocorrelation function provides us with valuable information 
on the dynamics of the trajectories of the system. It represents an important quantity 
strongly related to the physics described by the dynamical system. The construction of a 
dynamical system with prescribed statistical properties can have the background of either 
some experimentally given data for which a dynamical system should be modelled or it can 
be justified by the fact that a definite behaviour of the correlation function is required in 
the context of the control of a physical system ||^, [l^ . 



Restricting ourselves to the first part of the problem, i.e. the construction of a map with 
a given invariant density (the so called inverse Frobenius- Perron problem (IFPP)) we have 
recently found a general and for practical purposes very helpful representation of the solution 
of the IFPP within the class of smooth complete and unimodal maps 0. However the 
combined problem (both the invariant density and time autocorrelation function are given) 
is much more complicated. Some progress in this direction can be made if one uses the class 
of piecewise linear Markov maps as a basis for this search. Constructing such a map with a 
given autocorrelation function one can then obtain the proper invariant density performing 
a suitable conjugation (coordinate) transformation 0. A substantial disadvantage of this 
method is that the changes of the autocorrelation function in the course of the conjugation 
transformation are not fully under control. Another rather anaesthetic aspect is the fact that 
in many cases the solution found to the inverse problem does not fulfill certain smoothness 
criteria. 

The purpose of this letter is to present a new approach to the inverse problem including 
the autocorrelation function. It is based on the general representation of solutions of the 
IFPP found in ^ for complete smooth unimodal maps. We suitably parametrize the key 
function hj occurring in this approach and develop a stochastic (Monte-Carlo) algorithm to 
determine the optimal solution such that the x^ deviation of the resulting autocorrelation 
from a desired function is minimized. The paper is organized as follows: In section 2 we 
briefly review the representation of our solution to the IFPP. In section 3 we present the 
Monte-Carlo algorithm used to perform the x^ minimization with respect to the time auto- 



correlation function. Finally in section 4 we provide several applications in the framework 
of our Monte-Carlo scheme and discuss the convergence properties of our approach. 

2 The general solution of the IFPP for smooth uni- 
modal 1-d maps 

In a recent paper we investigated the problem of designing dynamical systems which 
possess an arbitrary but fixed invariant density. The solution to this inverse problem is 
of great interest for the numerical simulation of real physical systems as well as for the 
understanding of the relationship between the functional form of the map and the statistical 
features of the corresponding dynamics. We were able to derive a general representation 
of all ergodic and chaotic complete unimodal maps with a given invariant density. This 
corresponds to the general solution of the so called inverse Frobenius-Perron problem p|, |^, 
^ |T3|, |T^, [TB|, |T^ within this class of maps. We review here the basic aspects of our approach. 
As a starting-equation for the construction of the map we use the Frobenius-Perron 
equation: 

p{y)\dy\= Y. p{xi)\dxi\ (1) 

where the summation runs over all preimages of y (for unimodal maps i = L(left), R{right)). 
For any given complete unimodal f{x) the right preimage of y is determined if the left preim- 
age is given and vice versa. The essential feature of our approach is the following: A pre- 
scribed relation between the two preimages reduces the number of independent differentials 
on the rhs of eq.(P and allows to integrate the Frobenius-Perron equation. Such a relation 
is given by the function hf{x) which maps the position of the left preimage onto the position 
of the right one: 

xr = hf{xL) with /(xl) = /(x/j) (2) 

where Xmax is the position of the maximum of the map. hf{x) is a monotonously decreasing 
function on the defining interval and it is differentiable with the exception of a finite number 
of points. It obeys the equations: 



h'f{x) < Xe [0,Xmax] 

/i/(0) = 1 /i/(2;max) = 2;max (3) 



In terms of hf we obtain the general solution to IFPP within the class of unimodal maps as: 

fix) = fi-'{l-\fiix)-fiiHfix))\) (4) 

where Hf{x) is: 

^/w=(t-v\ ' '-^:'^- (5) 

[ hf [x) ; Xmax < a^ < 1 

rx 

and //(x) = / p{y)dy is the corresponding invariant measure. All unimodal maps with 
prescribed invariant density p{x) are given by (§), /i/(x) taking on all possible functional 
forms obeying (0) and (^ (for more details on the above derivation as well as its relation to 
the standard conjugation procedure we refer the reader to 0). 

The above shows that fixing the invariant measure is a relatively weak constraint in the 
framework of the inverse problem and that there is still a considerable freedom to model 
the mapping. In the next section we will use this freedom and suitably parametrize hf{x) 
such that, by tuning appropriately the corresponding parameters, we get a map with an 
autocorrelation function possessing a minimum deviation (to be defined below) from a desired 
correlation function. 

3 Stochastic optimization as a tool for the construction 
of a map 

Before presenting our search method for the optimal map in the sense of some desired 
statistical properties let us specify in more detail what we want to achieve and what our 
input for the problem at hand is. Consider a given invariant density p{x) (corresponding 
to a finite measure) which arises in the asymptotic limit (t -^ oo) from the dynamics of 
some unknown 1 — d fully chaotic single humped map. Consider also as given the first m 
values {C(l), C(2), .., C(m)} of the time autocorrelation function C{n) (note that C(0) is 
determined entirely through p{x)). We are seeking a map f{x) which possesses the above sta- 
tistical properties (p(x), C{n){n = 0, ..., m)). We focus in the following on the representation 
(^ for all admissible maps where the auxiliary functions hf{x) fulfill the requirements (§,|]). 
Using this expression for f{x) it is guaranteed that the map we are looking for possesses 
the invariant density p{x) (provided of course that p.{x) = Jq pit) dt in (§) exists). All our 
freedom in modelling the map f{x) is now contained in the auxiliary function hf{x). We will 
use a suitable representation of hf{x) to parametrize the map f{x). This parametrization is 
the starting point for a stochastic optimization procedure to obtain a map with the desired 
autocorrelation function C{n) at times n = 1, 2, .., m. One of the simplest ansatz would be 



to write down a piecewise linear expression for hf. Our knowledge of the properties and 
dynamics oi 1 — d maps however suggests that the local behaviour of the map in the neigh- 
bourhood of special points, like for example a marginal unstable fixed point in the case of 



intermittent dynamics (see ref . ^^ and references therein) , plays a crucial role in determining 
the statistical properties of the map. A more intuitive ansatz is therefore needed. Here we 
will proceed as follows. 

We define a 1 — d lattice with A^ + 1 points in the interval [0, Xmax]- The coordinates of the 
lattice points are Xp(z), z = 0, 1, .., A^ with Xp{0) = and Xp{N) = Xmax- The function hf{x) 
is piecewice defined. We will call in the following the expression of hf in the i — th interval 
[xp{i — l),Xp{i)) an element and use the notation hf^i for the i — th element. In practice 
one is free to choose many different expressions for the z-th element. Constraints on hf as 
for example continuity in [0, x^ax] together with the conditions (|^,0) however restrict the 
possible forms. A desirable but not necessary requirement is that hf^i should be analytically 
invertible in order to easily extend the solution in [0,a;max] to the interval (xmax, 1] (see (H)). 
Here we will investigate two different choices for the element hf^i. In both cases continuity 
of hf is guaranteed. The first case, refered to in the following as model I, is obtained by 
choosing the order of hf^i at the left point Xp{i — 1) in the interval [xp{i — l),Xp{i)). This 
parametrization allows us to fit the derivative of the map at a; = and has the advantage 
of using a minimal set of parameters in describing the element hf^i. The disadvantage of 
this choice is that one cannot fit independently the order of the maximum of the map at the 
right point Xp{N) of the last interval [xp{N — l),Xp{N)]. This is the reason for considering 
also a second case (model II). Here we use for the element hf^i an expression representing 
expansions around both limiting points Xp{i — 1), Xp{i) which match together at some point 
of the interval [xp{i — l),Xp{i)]. An additional parameter is therefore required in order to 
satisfy continuity for hf^i and its first derivative in {xp{i — l),Xp{i)). Due to this fact model 
II uses 2 more parameters than model I for the description of a single element /i/^j. 

Let us now discuss the two models in more detail. For model I the element hf^i is given 
by: 

hf,i = (?/p(0 - yp{i - 

X G [Xp{i — l),Xp{i)) ; i = l,2,..,N (6) 

where yp{i — l), Vpii) are the values of hf at the points Xp{i — 1), Xp{i) respectively, while the 
power a{i) determines the local behaviour of hf around Xp{i — 1). Due to the constraints 
(HJ^) we get: 2/p(0) = 1 and yp{N) = Xmax- The monotony of hf implies that yp{i — l) > Vpii) 
for i = 1,2,..,N. The above ansatz determines the map f{x) in the interval [0,a;niax]- To 
find f{x) in the remaining interval [xmax, 1] we take advantage of eq.(^) and therefore have 




to invert the auxiliary function hf{x). It is straightforward to show that the expression (^ 
can be inverted analytically leading to a closed form for hj^{x). The above defined ansatz 
for hf is in general piecewise smooth, i.e. smooth with the exception of a finite number of 
points. 

In model II we use the following ansatz for the element hf/. 

hf,i{x) = yp{i-l) -CL{i){x-Xp{i-l))°''''-'^ ; X e[Xp{i-l),Xs{i)) (7) 

hf,i{x) = yp{i) + CR{i){xp{i) - x)"""^'^ ; x e [xsii),Xp{i)] (8) 

where we have introduced four new parameters CL{i),CR{i),aji{i),Xs{i) and a new point 
Xs{i) G [xp{i — l),Xp{i)) for each i = I, ..., A^. Two of these parameters, namely cl(z) and 
CR{i), can be fixed demanding continuity of hf and its first derivative at x = Xs{i)- We arrive 
then at the following generating formula for hf-. 

hf^x) = yp{t - 1)- S, an(t){xsii) - Xp{i - l))i-"^«(x - Xp{t - 1))"^« ; 

X e [xp{i-l),Xs{i)) 

hf,i{x) = yp{i)+ S, aiW(2:pW-x,W)i-"««(Xp(2)-x)"«» ; 

X e [xs{i),Xp{i)) (9) 

with: 

^. (j/p(z - 1) - yp(z)) 

aR{i){Xs{i) - Xp{i - 1)) + aL{i){xp{i) - Xs{i)) 
hf{x) given in eq.(^ can also be inverted analytically to obtain h'f^ defined in [xmaxil]- 
Comparing eqs.(PJ^) we see immediately that the element of model II has 2 more parameters 
than the element of model I. 

In the remaining part of this section we will discuss the algorithm which allows us to 
obtain dynamical systems possessing certain statistical properties based on the above ansatz 
for the auxiliary function hf{x). To this end we will concentrate on the parametrization of 
model I. One can then directly apply these ideas to the case of model II. We have determined 
the ansatz of hf{x) and therefore also of f{x) in terms of the parameters {xp{i), yp{i), a{i)}. 
Once the values of these parameters are chosen the map is completely specified. Due to its 
appearance (see eq.(^)) we automatically know its invariant density (measure) and expecta- 
tion values of observables. The corresponding autocorrelation function can be obtained via 
its defining formula 

C{n) = I xf''\x)dix{x) - U xdnix)) (11) 

which requires a numerical integration. Here f^"'^ is the n — th iterate of the map /. 

To proceed with our central subject, the inverse problem, we assume that the first m 
values of the autocorrelation function are given. This can be due to some experimentally 



given data for which a dynamical system should be modelled or due to the fact that a 
definite behaviour of the correlation function is required in the context of the control of a 
physical system. As discussed above each set {xp{i),yp{i),a{i)}, i = 1,2,..,N respecting 
the constraint of monotony determines uniquely the autocorrelation function Chj {n) of the 
corresponding map /(x). One can now ask for the best set {xp{i),yp{i),a{i),i = 1,...,A^} 
in the sense that the resulting ChAn) possesses the least possible deviation from the given 
autocorrelation C{n). As a measure for the above-mentioned deviation one can use a x^-like 
cost function. We therefore introduce the following quantity: 



K[hf] 



tS(^^4^) 



The functional K[hf] is a highly nonlinear function of the parameters {xp{i),yp{i), «(i)} and 
we are looking for the global minimum of this function. To perform the minimization of 
K[hf] we use a Monte-Carlo (MC) approach based on the Metropolis algorithm |TH . 



The minimization is performed in several steps, increasing in each step the number of 
elements hf^i used for the determination of hj. We start with a lattice consisting of only two 
points (the origin Xp{0) = and the position of the maximum Xp{l) = Xmax)- This means 
that only one element hf^i is needed for the specification of hf. The parameters to be fitted 
in this case are only two: x^ax (3^p(l)) and the power a(l) determining the behaviour of 
hf{x) in the neighbourhood of the origin. The first step ends when the MC minimization 
has converged to some optimal values for the two fit-parameters. In the second step we use a 
lattice with three points 2:^(0), 2:^(1), Xp{2) with Xp{0) = 0, 2:^(1) = ^^, Xp{2) = Xmax- Now 
we need two elements /i/^i and hf^2 to determine hf. We do not keep the values of the old fit- 
parameters (xmax,tt(l)) fixed in the second step. Instead we use Gaussian distributed random 
variables for the choice of the old fit-parameters. The mean values of these Gaussians are the 
optimum values obtained for these parameters in the previous step and the corresponding 
widths are taken small enough to allow only weak fluctuations (~ 10%) around the mean 
values. Again we perform a MC optimization to obtain optimal values for the two new 
parameters. This procedure is repeated until the desired convergence is achieved. In each 
step the lattice size is increased by one point while we include two new parameters. 

We use the Metropolis algorithm to find the optimal values for Xp{i) and a{i + 1) in each 
step. For every trial in the i-th turn we assume that Xp{i) and a{i + 1) follow a uniform 
distribution in (0,1) and (0,oo), respectively (in practice the interval (0, oo) is replaced by 
the finite one (0,c) with an upper cutoff c >> 1). The annealing is introduced through a 
thermalized probability distribution of the type: P = e~^^^f^^'^ to avoid the trapping into 
local minima. The parameter T, playing the role of the temperature, is positive and has 



to be tuned adiabatically to smaller and smaller values such that the global minimum is 
reached asymptotically. 

We do not keep the values of the old fit-parameters fixed in the following step. Instead 
we relax this constraint using Gaussian distributed random variables for the choice of the 
old fit-parameters. The mean values of these Gaussians are the optimum values obtained for 
these parameters in the previous step and the corresponding widths are taken small enough 
to allow only weak fluctuations (~ 10%) around the mean values. 

Adding new fit parameters allows us to determine hf in more and more detail and im- 
proves the convergence of the autocorrelation function of the model dynamical system to 
the given (experimental) autocorrelation function. Indeed, as we shall see below, only a few 
elements hf^i are required to achieve a rather good convergence. In the next section we give 
some examples demonstrating how the above- described method can be applied to construct 
a (piecewise) smooth map simulating a system with given time correlations. 

4 Numerical examples and discussion 

Before turning to the examples and results of our computational method let us provide some 
additional aspects concerning the determination of the correlation function. We restrict our 
investigations to a rather small set of values for the correlation function C{n), more precisely 
to the set {C(0), ...,C(5)}. The reason for this restriction is that the exact calculation of 
the correlation function which has to be accomplished for each single step of the Monte- 
Carlo scheme is computationally very intensive (see below, the total amount of CPU of 
our calculations on a powerful workstation was approximately three months). The reliable 
evaluation of the correlation function is by no means trivial. As demonstrated for example in 
ref. |l^ the results obtained for the correlation function calculated with the trajectories of the 
dynamical system are, in many cases, not reliable and cannot be improved by going to longer 
propagation times. Therefore other methods for the calculation of the time correlations are 
needed. Here we use a numerical approach to eq.(0). It is based on the extraction of 
the monotony intervals for the n-th iterate. The endpoints of the monotony intervals are 
given as the preimages of the maximum Xmax- The integration is then performed for each 
monotony interval separately. This ensures an accurate although very CPU time consuming 
evaluation of the correlation function. The latter is related to the exponential growth (2") 
of the number of monotony intervals with increasing n. 

We apply now the stochastic method described in the previous section to several examples 
to demonstrate its capability of providing dynamical systems with prescribed statistical 
properties. The four different cases we study here can be summarized with respect to their 



statistical properties as follows: 

• exponentially decaying autocorrelation function and uniform density 

• exponentially decaying autocorrelation function and linear density 

• oscillatory decaying autocorrelation function and power-law density 

• power-law decaying autocorrelation function and uniform density 

Figs. (1-4) illustrate the corresponding results of our Monte-Carlo approach. The subfig- 
ures (a,c) within each figure show the resulting maps for model I and model II, respectively. 
Each of the subfigures (b,d) within each figure (1-4) shows both the prescribed data for the 
autocorrelation signal as well as the result of our optimization approach for model I and II, 
respectively. In comparing the given autocorrelation data and the results of our optimization 
scheme we deduce in the following a 'mean relative error' per individual data point l^ of 
the autocorrelation function. 

Since it is not our goal to provide as precise data as possible for the optimized autocorre- 
lation functions but to demonstrate the feasibility of our approach with a good accuracy for 
a few data points we restrict ourselves to auxiliary functions hf composed of a few elements 
/i/ j. For the case of model I we use (with one exception) four points, i.e. three elements 
hf^i for the decomposition of the auxiliary function hf. For model II we use only two points, 
i.e. one element hf^i. Typically a few ten thousand Monte-Carlo steps are performed. The 
most time consuming part of the algorithm is the calculation of the autocorrelation function 
which has to be done for each MC step. 

First we consider a simple example for which the maps we are looking for are well- 
known. We use an exponential decay as a prescribed behaviour for the autocorrelation 
signal while the invariant density of the dynamical process is assumed to be uniform in 
[0, 1]. A corresponding family of maps, the nonsymmetric tent maps, has been studied in 
some detail in ref.[^]. They are the nonsymmetric tent maps. The resulting maps of our 



Monte-Carlo optimization are illustrated for model I in fig. 1(a) and for model II in fig. 1(c). 
They show only minor differences, i.e. the outcome of approach I is almost the same as for 
approach II, and both are to a very good approximation an asymmetric tent map. Both 
models lead also to a very good approximation of the corresponding prescribed correlation 
function with an error of only 2.2%. 

For the second example the two models lead to quantitatively different results. The 
invariant density is supposed to be linear while the autocorrelation function shows an expo- 
nential decay similar to the previous case. The results of the stochastic minimization for this 
case are illustrated in figure 2. The two obtained dynamical systems (figs.2(a,c)) have a very 



different appearance reffecting the fact that our modelhng procedure allows for various dy- 
namical systems with the same invariant density and autocorrelation. Let us briefly address 
the main features of the resulting maps. The map of model I possesses an almost marginal 
unstable fixed point at the origin and an almost vertical derivative at a; = 1. Additionally 
it possesses two obvious points of noncontinuous derivatives located on the left and right 
branch of the map, respectively. The one on the left branch has a nonvanishing right deriva- 
tive which, via the function hf{x), is mapped to a point with almost vertical left derivative 
in the right branch of f{x). The reader should note that the latter point coincides with the 
nonzero fixed point of the map ! The map of model II (fig.2(c)) is almost straighlined on 
its left branch and shows an almost vertical right derivative at the maximum and an almost 
vertical left derivative at x = 1. In general the map of model II is 'smoother' than the map 
of model I, which is an overall tendency to be observed in any of our examples. It can be 
viewed as a result of the additional flexibility within model II which allows to independently 
adapt the left and right derivative in a given interval thereby joining them smoothly together 
(see the above description for the ansatz of hf in model II). The prescribed and optimized 
autocorrelation data for model I are illustrated in subfigure 2(b): They show a deviation of 
129% which is predominantly due to the inability of reproducing the single point C(l) within 
the approach of model I. Although hardly visible in fig. 2(b) the prescribed and optimized 
data coincide very well for {C{n), 2 < n < 5}. The optimized correlation function for model 
II leads to a much better approximation to the prescribed exponential decay and yields an 
error of only 5.4% (see fig. 2(d)). 

As a third example we use an invariant density obeying a power law with an exponent P = 
—0.5 while the autocorrelation function is chosen to possess a decaying oscillatory behaviour. 
The corresponding minimization results are shown in fig.3. Both maps (figs.3(a,c)) show an 
almost horizontal derivative at a; = 1 and a strong cusp at ajmax- For the map of model II 
the left derivative at the maximum is almost vertical. Again the appearance of the map of 
model II is much smoother compared to the map of model I. Regarding the autocorrelation 
function we have the opposite situation compared to the previous example. Within model 
I (fig.3(b)) we obtain a relative error of 22% for the autocorrelation data while model II 
(fig. 3(d) provides an error of 45%. Obviously model I is advantageous in the present case (at 
least for the present number of grid points chosen for the ansatz of hf): it nicely reproduces 
the oscillations of the correlation function. 

Finally we study the case that the autocorrelation function decays algebraically with the 
exponent 7 = —2.5, i.e. we encounter the case of long-range correlations. The invariant 
density is chosen to be uniform. Both resulting maps (fig.4(a,c)) possess a large but finite 
derivative at x = 1. The map of model I possesses an almost vertical derivative at a single 
point on its right branch. In contrast to our second example (see above) this point does not 
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coincide with the fixed point. It is conjectured that this single point with almost vertical 
derivative is responsible for the observed power law decay of the correlation function (fig4(b)) 
(see also ref.|]21[]). Analogously the map of model II (fig. 4(c)) shows a very well-pronounced 
vertical derivative at the cusp. It is interesting to observe that in this case both models 
I and II provide very satisfactory approximations to the prescribed autocorrelation data 
(fig4(b,d)), i.e. a relative error of only 13%. 

5 Summary 

We have introduced a stochastic Monte-Carlo based approach to the inverse problem which 
uses both the invariant density as well as a finite number of points of the autocorrelation 
function as prescribed statistical quantities. A key ingredient for this approach was a recently 
found representation for one- dimensional complete chaotic and single-humped dynamical 
system in terms of an auxiliary function hf. This representation is a formal solution of 
the inverse Frobenius- Perron problem and provides the dynamical system f{x) explicitly as 
a function of the given measure and the function hf. hf therefore refiects the freedom of 
changing the map without changing its invariant measure (density). In order to quantify the 
freedom available for the determination of hf, i.e. to parametrize its functional space, we 
have introduced two different models which allow to vary hf extensively. The parameters 
involved are then used within our stochastic minimization procedure to obtain a correlation 
function with least deviation from a prescribed autocorrelation signal. Through a number of 
examples we have demonstrated that our approach possesses an enormous fiexibility allowing 
for a large variety of qualitatively different behaviour of the density and correlation function. 
To our knowledge there is in general no unique map which belongs to a given density and 
correlation function. This fact has to be seen in the context of the present investigation as 
an advantage since it allows for a great flexibihty and possible variety with respect to the 
underlying dynamical systems. 

We would like to mention that the above-discussed prescribed behaviour of the auto- 
correlation data (exponential, oscillating, power law decay) is, strictly speaking, enforced 
only for the first five points included in our Monte-Carlo optimization. In principle it is 
imaginable that this behaviour represents a transient and the asymptotics of the correlation 
function might show a different behaviour. To determine exact asymptotic properties of cer- 
tain dynamical systems is however not the issue of the present paper. Our goal is to extend 
the inverse problem by including the correlation function in terms of a few (experimentally) 
available data points, thereby enabling us to design a dynamical system with desired statis- 
tical properties. Furthermore our approach might be suggestive in terms of influencing or 
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controlling dynamical systems [^ |T0| in a certain way 
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FIGURE CAPTIONS 

Figure 1: The stochastic minimization results for a dynamical system with uniform in- 
variant density and exponentially decaying autocorrelation function, (a) The resulting map 
using model I. (b) The autocorrelation function for the map of model I (solid line) and the 
corresponding prescribed data (full circles), (c) The resulting map using model II. (d) The 
autocorrelation function for the map of model II (solid line) and the corresponding prescribed 
data (full circles). 

Figure 2: Same as in figure 1 but for prescribed linear invariant density and exponentially 
decaying autocorrelation function. 

Figure 3: Same as in figure 1 but for prescribed power-law invariant density (/9 = —0.5) 
and oscillatory decaying autocorrelation function. 

Figure 4: Same as in figure 1 but for prescribed uniform invariant density and power-law 
decaying autocorrelation function (with exponent -2.5). 
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